Domain growth by isothermal aging in 3d Ising and Heisenberg spin glasses 
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Non-equilibrium dynamics of three dimensional model spin glasses - the Ising system 
Feo.soMno.soTiOa and the Heisenberg like system Ag(ll at% Mn) - has been investigated by mea- 
surements of the isothermal time decay of the low frequency ac-susceptibility after a quench from 
the paramagnetic to the spin glass phase. It is found that the relaxation data measured at different 
temperatures can be scaled according to predictions from the droplet scaling model, provided that 
critical fluctuations are accounted for in the analyzes. 
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A spin glass ages - the magnetic response of the sys- 
tem depends on its thermal history and especially the 
length of time it has been kept at a constant temper- 
ature in the glass phase Q]. An aging phenomenon is 
also common to many other disordered and slowly relax- 
ing physical systems 0|. A thorough understanding of 
the driving forces behind the aging process is thus im- 
portant to adequately model glassy matter. Aging in 
spin glasses (SG) is within real space SG models asso- 
ciated with a slow growth of the coherence length for 
equilibrium SG order. The droplet theory || provides 
concrete predictions of the scaling properties of the age 
dependent macroscopic dynamical susceptibility in terms 
of a growing coherence length Lr{t). One key test of the 
applicability of real space SG models is thus to quanti- 
tatively confront age dependent susceptibility data with 
the droplet model predictions. A positive result of such 
a study does not only give numbers for the parameters 
involved but allows further quantitative investigations of 
important ingredients of real space pictures such as chaos 
with temperature and overlap length, and could even pro- 
vide hints on the topological character of the assumed 
SG equilibrium domains. It does on the other hand not 
tell anything about the relevance of other approaches to 
spin glass physics. In this Letter, we investigate the time 
decay of the low frequency ac-susceptibility at constant 
temperature of two model spin glasses - the 3d Ising sys- 
tem Feo.5oMn .5oTi03 and the 3d Heiscnberg-like Ag(ll 
at% Mn) (AgMn). In contrast to the outcome of all 
earlier corresponding studies (see for example [§, ||), we 
here find that the relaxation data can be scaled according 
to the explicit relation between the growth of Lx(t) and 
the relaxation of the ac-susceptibility suggested by the 
droplet theory || , provided however that critical fluctu- 
ations are accounted for in the analyzes. 

In the droplet model it is assumed that isothermal ag- 
ing at a constant temperature T after a quench from a 
temperature above the spin-glass transition temperature 
T g is a coarsening process of domain walls as in many 
other phase ordering systems. During the isothermal ag- 
ing of a spin glass of spatial dimension d, the temporal 



ac-susceptibility at a given angular frequency lo at time 
t after the quench is supposed to scale as & , 
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where Lxit) is the typical size of the domains and 
Lt(1/<jj) is the typical size of the droplet being polarized 
by the ac field. The basic physical idea behind this scal- 
ing form is that the presence of a domain wall effectively 
reduces the excitation gap T(T)(L/Ln) e of droplet exci- 
tations that are touching the domain wall. Here 6 is the 
stiffness exponent and Lq is a microscopic length scale. 
The equilibrium ac susceptibility Xcq^) i s proposed to 
scale as JtJ, 
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where is an universal constant and q m is equal to the 
Edwards- Anderson order parameter for a symmetric spin 
glass. The dynamics is governed by thermally activated 
processes and the length scale L associated with a certain 
temperature T and the time scale t is given by 
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with tq(T) being the unit time scale for activated pro- 
cesses. This growth law is derived from scaling of the 
typical free-energy barrier associated with the length 
scale L, F barricr (L) ~ A(T)(L/L )^ with the exponent 
ip > 0. The stiffness for the excitation gap T(T) and 
barrier A(T) depends on temperature as T(T) ~ Je 8 " 
and A(T) ~ Je^ v , respectively, with J being the charac- 
teristic energy scale and e = T/T g — 1. 

An important point is that the unit time scale tq(T) 
for activated hopping processes is not simply given by 
the real microscopic time scale ||, which is r m ~ h/J ~ 
10~ 13 s in spin systems. A plausible choice for tq(T) is 
instead the critical correlation time r c as proposed in || , 
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with z and v being the dynamical critical exponent and 
the exponent for the divergence of the correlation length 
£(T), respectively. 

Due to a possible slow crossover from critical dynam- 
ics at t <C t c (T) to activated dynamics at t 3> t c (T), the 
logarithmic domain growth law Eq. (||) may apply only 
for ideally asymptotic regimes beyond the time scales of 
numerical simulations j9). Results from numerical sim- 
ulations [^0) accordingly suggest an alternative growth 
law for the spin glass coherence length 



L T (t) ~ L (i/r m ) 1/z(T) 



(5) 



where z(T) is a temperature-dependent effective expo- 
nent. Inserting Eq. (^) into the basic scaling forms 
Eq. ([l]) and Eq. (§) one can obtain corresponding for- 
mulas for the relaxation and frequency dependence of x" 



X "{LO,t)- X '^) 
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with temperature dependent effective exponents 



(6) 



fc(T) = (d-0)/z(T), a(T) = 9/z(T). (7) 

Interestingly, Eq. (g) is very close to that of the cut - 
scaling proposed in previous experiments j5j where in- 
deed the exponents are found to depend on temperature. 
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FIG. 1: x" (^) vs time for the Feo.soMno.soTiOa sample 
with lu/2it = 0.17,0.51,1.7,5.1,17,55,170 Hz (from bottom 
to top). The left inset show x"(w) vs temperature, measured 
on cooling, for uj/2tt = 0.17, 5.1, and 170 Hz (from bottom to 
top). The right inset show x"^) vs time for T = 15, 17, and 
19 K (from bottom to top). 

The experiments were performed on a single crystal 
of the short-range Feo.5oMn .5oTi03 Ising spin glass Q 
and a polycrystalline sample of the long-range AgMn 
Heisenberg-like spin glass. The isothermal aging was 



measured by recording the ac-susceptibility, as a func- 
tion of time (for up to 10 4 s), after rapidly cooling the 
spin glass from a temperature above T g to the measure- 
ment temperature. The amplitude of the probing ac field 
was 0.02 Oe and frequencies in the range uj/2ir — 0.17 - 
170 Hz were used. The background field was smaller than 
0.1 mOe. The cooling rate, 0.05-0.08 K/s, was the fastest 
allowed by the non-commercial squid magnetometer E^] 
employed in these studies. Some typical experimental re- 
sults for the Feo.soMno.soTiOa sample are shown in Fig. 
[I]. The relaxation toward equilibrium is for both samples 
so slow, even for the highest measurement frequencies, 
that Xeq( w ) could not be estimated directly but had to 
be included as a fitting parameter in the analyzes at all 
frequencies and temperatures. Due to the finite cooling 
rate the system will have an "effective age" when it has 
been cooled to the measurement temperature. An "effec- 
tive time" t e ff has therefore been added to the measure- 
ment time in the analyzes. In the following analyzes we 
have used t e g(T = 30, 25, 20 K) = 10, 15, and 30 s for the 
AgMn sample and t cS {T = 19, 17, and 15 K) = 150, 200, 
and 350 s for the Feo.5oMn .5oTi03 sample M. 

The isothermal aging data were first analyzed assum- 
ing the logarithmic domain growth law in Eq. (0) , which 
inserted in Eq. ([j]) yields, 
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This scaling law should permit all isothermal aging data 
from measurements at different temperatures and fre- 
quencies to collapse onto one single master curve. A 
necessary condition for such a scaling is activated dy- 
namics, i.e. t > (1/w) 3> t c {T). This condition is not 
fulfilled at T = 19 K for the Feo.5oMn . 50 Ti0 3 sample, 
where t c = 9 • 10~ 4 s which is equal to l/u> for the highest 
frequency. The data at T = 19 K have therefore been ex- 
cluded from this analysis. The parameters describing the 
critical dynamics were varied within reasonable empirical 
values (T g = 32.5 - 33 K, zv = 7 - 8, and r m w 10" 13 s, 
for the AgMn sample; T g = 21.0-21.5 K, zv = 10-11, 
and r m rj 10~ 13 s for the Feo.5oMn .5oTi03 sample). For 
the AgMn sample, a good scaling could be obtained us- 
ing critical parameters in the indicated range resulting in 
(<f-0)/V> w 1.7 [c.f. Fig.|(b)]. For the Fe .5oMno. 5 oTi0 3 
sample the quality of the scaling was more sensitive to the 
parameters determining the critical dynamics. However, 
a reasonably good scaling was obtained and we can esti- 
mate that (d — 9)/ip ~ 1.5 [c.f. Fig. ||(a)]. It was not pos- 
sible to scale the data for any of the two spin glasses using 
a constant value for To ~ 10 -13 s. The remaining fitting 
parameter in the scaling analysis, x'cqi^^T), is predicted 
to show a frequency and temperature dependence accord- 
ingtoEq. (§), i.e. cx [Tln((l/^)/r c (T))]-( 1 +^). 

Asymptotic fits of the obtained Xg q (u;,T) toward the 
highest values of Tln[(X/u))/r c (T)] to this equation then 
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ln[t/x c (T)]/ln[(1/(o)/T c (T)] 

FIG. 2: [x"M) - X'M]/X"(u,t) vs 

ln[t/r c (T)]/ln[(l/cj)/r c (T)] on a log-log scale (a) 
for the Feo.soMno.soTiOs sample with u>/2ir = 
0.17,0.51,1.7,5.1,17,55,170 Hz, T = 15,17 K and (b) 
for the AgMn sample with lo/2tt = 0.17,1.7,17,170 Hz, 
T = 20,25,30 K. The insets show x" q M vs ln[(l/u)/r c (T)] 
on a log- log scale (a) at T = 17 K (pluses), and 15 K (circles) 
and (b) T = 30 K (triangles), 25 K (crosses), and 20 K 
(squares). 



give 1+0 ftp w 1.8 for the AgMn sample and 1+0 /ip ~ 1.1 
for the Feo.5oMn .5oTi03 sample. Combining the results 
for the exponents we obtain ip rj 1.2 and w 1.0 for 
the AgMn sample while ip m 1.9 and 6* w 0.2 for the 
Feo.5oMn .5oTi03 sample. These estimates are within 
the imposed limits < (d — l)/2 and < ip < d — 1 
j/J . For both samples, the exponent (d — 0) /ip is strongly 
dependent on the parameters describing the critical dy- 
namics and increases with decreasing t c (T) (i.e. higher 
T g and/or smaller zv), while the exponent 1 + /ip was 
less sensitive. The here obtained value of ip for the AgMn 
is close to a value reported earlier tp « 1.3 p4| , and the 
extracted value of for the Ising spin glass is in accord 
with estimations from numerical simulations on the 3d 
Edwards- Anderson model j|, |l5| . On the other hand, for 
the Ising spin glass ip rs 1.9 is much larger than earlier re- 



ported experimental values ip « 0.8 jl6| and ip w 0.3 — 0.7 

To give some further perspective on the values ob- 
tained for the exponents and ip, it should be pointed 
out that employing the alternative form [x"{uj,t) — 
Xeq( w )]/Xeq( w ) 011 the right hand side of Eq. (0) yields an 
almost equally good scaling behavior, but that the val- 
ues of the exponents alters to ip R* 1.3 and 0.13 for 
the Feo.5oMn .5oTi03 system and ip rs 0.8 and w 0.6 
for the AgMn sample. Here it is worth to note that the 
alternative scaling form and the scaling form of Eq. (|l|) 
are expected to give the same result if correction terms at 
order (Lt(1/<jj)/ LT(t)) 2 ( d ~ s ' are negligible. However the 
ratio Lt{1/^j)/ Lrit) cannot reach values much smaller 
than 1 on experimental time scales so that neglecting this 
term might have affected the analyzes. Also, the possible 
slow crossover from critical to activated dynamics may 
require a more complete functional form than only the 
asymptotic form Eq. (||) to precisely describe the growth 
law on experimental time scales. 

Assuming instead that the domain growth is cor- 
rectly described by the algebraic growth law Eq. (0), 
the isothermal aging should scale according to the power 
law in Eq. (||) with exponents defined in Eq. (Q). In 
this case, the data corresponding to different frequencies 
could be collapsed into one master curve for each tem- 
perature (see Fig. |J), using Xg ? (w,T) as fitting param- 
eter. The Xeq^) values obtained are shown in log- log 
plots in the insets of the figure together with fits to the 
power law behavior predicted in Eq. (^]). The tempera- 
ture dependence of the exponent k{T) [c.f. Eq. (R)] ac- 
cords roughly with the behavior suggested by numerical 
simulations z(T) ~ z(TjT) 0. However, some ambi- 
guity appears when examining the ratios a(T)/k(T) of 
the derived exponents (see inset of Fig. |^). At least for 
AgMn this ratio clearly increases with decreasing tem- 
perature, while Eq. (Q) implies that it should be constant 
and equal to 0/(d — 0). In the case of Feo.5oMn .5oTiC>3 
the derived rj 0.65, whereas for AgMn appears to 
increase from about 0.9 at 30 K to 1.2 at 20 K, which is 
outside the imposed limit < 1. 

It has, as mentioned in the introduction, earlier been 
found that non-equilibrium relaxation in spin glasses can- 
not be described by ln(t)-scaling (see for example Refs. 

[5] ) . However, in those investigations it was always as- 
sumed that tq(T) = T m . Quite remarkably, we now find 
that, accounting for critical fluctuations by identifying 
To(T) with the critical correlation time, the ln(t)-scaling 
permits to collapse all data measured at different fre- 
quencies and temperatures onto one single master curve 
(see Fig. ^) . The importance of including critical fluctu- 
ations when analyzing non-equilibrium dynamics in spin 
glasses has also been emphasized in recent experimental 
work by E. Vincent and co-workers ]l4| |. It can further 
be mentioned that numerical studies of the 4d Edwards- 
Anderson model only allow a logarithmic growth law if 
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FIG. 3: tx"(w,t) - >&(w)]/x"(w.*) vs wt on a log-log 
scale (a) for the Feo.soMno.soTiOs sample with u/2n — 
0.17, 0.51, 1.7, 5.1, 17, 55, 170 Hz and (b) for the AgMn sam- 
ple with lo/2-it = 0.17, 1.7, 17, 170 Hz. The left inset in (a) 
shows a(T)/k(T) vs T/T g for the AgMn sample (circles) with 
T g — 32.8 K and for the Feo.soMno.soTiOs sample (triangles) 
with T g — 21.3 K. The other insets show Xeq^) vs w on a 
log- log scale (a) at T = 19 K (stars), T = 17 K (pluses), and 
15 K (circles) and (b) T = 30 K (triangles), 25 K (crosses), 
and 20 K (squares). 



critical fluctuations are accounted for [jfj. 

The current finding, that the values of the exponents 9 
and ij) for the Feo.soMno.soTiOs Ising spin glass are sig- 
nificantly different from those of the AgMn Heisenberg- 
like spin glass, is consistent with the chiral glass picture 
proposed by Kawamura JT7fl , conjecturing different uni- 
versality classes for Heisenberg and Ising spin glasses, 
and numerical simulations within that picture indicate, 
in accord with our results, that the value of should be 
larger for a vector spin glass than for an Ising system 
jig . Further experimental support for the existence of 
different universality classes in spin glasses was recently 
reported from magnetic torque measurements, yielding 
significantly different critical exponents for isotropic and 
anisotropic spin glasses JIIJ. 

The numerical values of exponents involved in the 



growth law of the coherence length that have been ob- 
tained in this study will allow quantitative investigations 
of other aspects of spin glass dynamics within the droplet 
model, such as analyses of chaos with temperature and 
overlap length p(i[ . 

We thank Eric Vincent and Roland Mathieu for useful 
discussions. This work was financially supported by the 
Swedish Research Council (VR). 
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